Can there be neutrino oscillation in Gamma-Ray Bursts fireball ? 
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The central engine which powers the Gamma-Ray Burst (GRB) fireball, produces neutrinos in the 
energy range of about 5-20 MeV. Fractions of these neutrinos may propagate through the fireball 
which is far away from the central engine. We have studied the propagation of these neutrinos 
through the fireball which is contaminated by baryons and have shown that, resonant conversion of 
neutrinos are possible for the oscillations of v e <-» v^r, v e *-* v s and z/ MiT <-+ v s if the neutrino mass 
square difference and mixing angle are in the atmospheric and/or LSND range. On the other hand 
it is probably difficult for neutrinos to have resonant oscillation if the neutrino parameters are in 
the solar neutrino range. From the resonance condition we have estimated the fireball temperature 
and the baryon load in it. 

PACS numbers: 14.60.Pq,98.70.Rz 



GRBs are short, non-thermal bursts of low energy (~ 
100 KeV-1 MeV) photons and release about 10 51 -10 53 
ergs in a few seconds making them the most luminous ob- 
ject in the universe. A class of models call fireball model 
seems to explain the temporal structure of the bursts 
and the non-thermal nature of their spectraQ, 0, 0, S @ 
Sudden release of copious amount of 7 rays into a com- 
pact region with a size cSt ~ 100-1000 KmpjJ creates an 
opaque 7 — e fireball due to the process 7+7 ~* e + + e _ . 
The average optical depth of this process[6| is r 77 ~ 10 13 . 
This optical depth is very large and even if there are no 
pairs to begin with, they will form very rapidly and will 
Compton scatter lower energy photons. Because of the 
huge optical depth, photons can not escape freely. In the 
fireball the 7 and pairs will thermalize with a temper- 
ature of about 3-10 MeV. The fireball expands relativis- 
tically with a Lorentz factor V ~ 100 — 1000 under its 
own pressure and cools adiabatically due to the expan- 
sion. The radiation emerges freely to the inter galactic 
medium (ISM), when the optical depth is t 77 ~ 1. 

In addition to 7, e pairs, fireballs may also contain 
some baryons, both from the progenitor and the sur- 
rounding medium. These baryons can be either free or 
in the form of nuclei. If the fireball temperature is high 
enough (more than 0.7 MeV), then it will be mostly in 
the form of neutrons and protons. Derishev et al. 0, 13 
argued that roughly equal numbers of neutrons and pro- 
tons should be present in the fireball. The electrons asso- 
ciated with the matter (baryons) can increase the opac- 
ity, hence delaying the process of emission of radiation 
and the baryons can be accelerated along with the fire- 
ball and convert part of the radiation energy into bulk 
kinetic energy. But irrespective of it, the baryonic load 
has to be very small, otherwise, the expansion of the 
fireball will be Newtonian, which is inconsistent with the 
present observations. Why the baryonic loading is so low 
in the fireball is still an open question to be answeredQ. 
The neutrino oscillation may overcome the baryon load- 



ing problemjij E3 ' lne evolution ot t he p ure 
(with no baryons) has been studied in ref. [ll], [T^] • In the 
expanding fireball, protons are accelerated in shocks and 
collide with the photons to produce charged pions, which 
give rise to ultra high-energy neutrinos [13j. These neu- 
trinos can be detected by km 2 detectors. Observation of 
these neutrinos can study oscillating neutrino flavors in 
the largest possible baseline, test the equivalence princi- 
ple and many other neutrino properties fl4| . 

The hidden central engine which powers the fireball 
is still unknown, but observation suggests that it must 
be compact. The prime candidates are merger of neu- 
tron star with neutron star (NS-NS), black hole- neutron 
star binaries (BH-NS), hypernova/collapsar models in- 
volving a massive stellar progenitor Q, a^yiuJ The 
recent observations of GRB030329 strongly favor collap- 
sar model 16]. In all these models, the gravitational 
energy is released mostly in the form of uv, gravitational 
radiation and a small fraction (~ 10 -3 ) is responsible 
to power GRB. Neutrinos of energy about 5-20 MeV are 
generated due to the stellar collapse or merger event that 
trigger the burst. Also due to nucleonic bremsstrahlung 
NN — » NNvP and e + e~ — * vv processes muon and tau 
types neutrinos can be produced|l7| during the merger 
process and its flux will be very small. Fractions of these 
neutrinos may propagate through the fireball which is far 
away from the central engine. Also within the fireball, 
because of the weak interaction process p + e~ — > n + v e , 
MeV neutrinos can be generated and propagate through 
it. The fireball plasma being in an extreme condition, 
may affect the propagation of MeV neutrinos through it. 

In a heat bath the dispersion relation which governs 
the propagation of the particle gets modified and this 
can have a drastic effect on the particle propagation in 
the heat bath. The neutrino pr opagati on in a thermal 
bath has been studied extensively 1 181 Il9l 120. 21] . In gen- 
eral, the propagating neutrino will experience an effective 
potential due to the particles in the heat bath. Particu- 
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larly the neutrino propagation in the early universe hot 
plasma, as well as in the supernova medium has profound 
implications in the respective physics. For example the 
neutrino oscillation in the early universe hot plasma may 
change the relative abundances of v e and v e and affect 
the primordial nucleosynthesis of light elements. In the 
supernova case, the neutrino oscillation in the dense and 
compact environment can affect the cooling mechanism. 

The electron type neutrinos have charge current as 
well as neutral current interactions, but the muon and 
tau types will experience only the neutral current in- 
teractions. To the leading order the effective potential 
experience by the neutrinos is proportional to the differ- 
ence of the particle anti-particle number densities. So 
if the system under consideration has equal number of 
particles and anti-particles, the leading order contribu- 
tion vanishes. On the other hand, the next-to-leading 
order contribution (i.e. the term proportional to 1/M 4 , 
where M is the vector boson mass) is proportional to the 
sum of the particle and anti-particle number densities 
and this is the leading contribution to the effective po- 
tential. The early universe hot plasma, supposed to have 
almost equal population of particles and anti-particles, 
hence the leading contribution to the neutrino effective 
potential will be proportional to 1/M 4 . Similar situa- 
tion can arise for neutrinos propagating in the GRB fire- 
ball if one consider mostly photon- lepton fireball|l|, 0) - 
Here in the present work we consider a fireball containing 
mostly photon-lepton with little baryon contamination, 
which mimics the early universe hot plasma and study 
the neutrino propagation within it. 

In a relativistic and nondegenerate e , proton and neu- 
tron plasma, the effective potential experiences by v e is 
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and by a muon or tau neutrino it is given by 
where 
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The asymmetry of the particle a is defined as 

(N a - N a ) 



L a = 
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where 7V 7 = -^^(3)T 3 is the number density of pho- 
tons. For the anti-neutrinos, the effective potential will 



be given by changing L a — > — L a . For a pure 7 and 
fireball, we have L p — L n = 0. 

We consider the neutrino oscillation processes: v e <-> 
v e *-> v s and the antineutrino processes. The ef- 



fective potential difference for v e 
processes, is given by 



<-> and v e v 



V ~ 4.02 x lQ- 12 Tl IeV [±L e ~ 6.14 x l^T 2 MeV ] MeV, 

(6) 

where ± corresponds to v and v respectively, and for 
v e <-* v s oscillation, the effective potential difference is 
given by Eq.JQ). For pure 7, e fireball (CP-symmetric), 
L e — and only the higher order term will contribute. 
The conversion probability for the above processes for a 
constant V is given by, 



, A 2 sin 2 29 . 7 ,uot s 
= 5 "in M. 



where 



(7) 



= \J(V- Acos26») 2 + A 2 sin 2 26», (8) 



with A = Am 2 /2i?„, V is the potential difference, E v is 
the neutrino energy and 9 is the neutrino mixing angle. 
The oscillation length is given by 
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where L v — 2ir/A is the vacuum oscillation length which 
can be recovered for V — 0. For resonance to occur we 
should have from Eq.JSJ) 



V = A cos 20. 



(10) 



We have so far assumed that lepton asymmetry in the 
fireball does not vary with distance. But in reality, the 
lepton asymmetry changes with distance. So in this case 
we will consider the adiabaticity of the resonant conver- 
sion. The adiabaticity condition at the resonance |22j is 
given by 



K r = 2.0 x 10~ 3 



At 
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dx 



> 1, 
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where l cm is some length scale expressed in centimeter, 
x is a dimensionless variable, Am 2 is in unit of eV 2 and 
TmeV and EmcV are expressed in unit of MeV. The above 
condition depends on the neutrino parameters and the 
change in the L e as one goes away from the center of the 
fireball. 

As already stated in the introduction, baryon loading 
is still an outstanding problem. It is believed that the 
contamination is very low (10 _8 A/ Q - 10~ 5 M Q )[l|,|2, so 
that, the fireball can have a ultra-relativistic expansion. 
For simplicity, we consider a charge neutral spherical fire- 
ball (L e = L p ) of initial radius R with equal number of 
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protons and neutrons in it. Then the baryon load in the 
fireball is 



M h 



2.23 x 10 



^rT MeV L e MQ, 



(12) 



where R 7 is in unit of 10 7 cm. As stated above, the 
baryon contamination is of order 1O _8 M0 to 1O _5 M0, 
which corresponds to L e ~ 4.47 x 10~ 5 i?f 3 T^ 3 y to 
L e ~ 4.47 x 10~ 2 i?f 3 T M 3 eV respectively. The thermal- 
ized fireball has a temperature T ~ (L/inaR 2 ) 1 / 4 ~ 
3-10 MeV, with L the luminosity and a the Stephen- 
Boltzmann constant. 

The processes (active to active oscillations) do not de- 
pend on the baryon asymmetry in the fireball, simply 
because, the neutral current contribution to the poten- 
tial of all the active neutrinos is the same and for active 
to active oscillation this contribution cancels out, leav- 
ing only the dependence on L e . But the active-sterile 
oscillation depends on both the lepton and the baryon 
asymmetry. In Eq.lJBj, V > for neutrino process if 



L E > 6.14 x 10- g T^ eV , 



(13) 



and depending on the fireball parameters and the neu- 
trino properties the resonance condition can be satisfied. 
On the other hand, antineutrino process will never sat- 
isfy the resonance condition because the potential is al- 
ways negative. The resonance condition of Ea. (|10l) can 
be written as 



L e T 3 MeV = 0.124 



Am 2 cos 29 



EMeV 

and the resonance length is given by 

T O/l Q E MeV 

L res ~ 248 — = cm. 

Am 2 sin 29 



(14) 
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Putting the value of L e from Eq. I|14|) in Eq. I|13|) , the con- 
straint on the fireball temperature is 



T MeV < 0.2 x 10 8 



Am 2 cos 29 



E 



MeV 



(16) 



Thus the fireball temperature derived in Ea. l|16|l is the 
one required in order to have resonant conversion of the 
neutrinos. As we have already discussed, during the stel- 
lar collapse or merger events that trigger the burst, neu- 
trinos of energy about 5-20 MeV are copiously produced 
and some of these neutrinos will propagate through the 
fireball. Apart from this, due to inverse beta decay, MeV 
neutrinos can also be generated within the fireball. So 
here we will take E-^eV = 5 and 20 for our calculation 
to estimate the fireball parameters. 

Let us study the resonance condition for Solar, At- 
mospheric and the reactor (LSND) neutrinos, where we 
know approximately the neutrino mass square differences 
and the neutrino mixing angles from the recent experi- 
mental results. These can constraint the fireball param- 
eters. 



The recent analysis of the salt 
combines with the KamLAND 



hase data of SNO|23j 
reactor antineutrino 



results gives 6 x 10 eV < Am 2 < 10 eV and 
0.8 < sin26> < 0.98 with a confidence level of 99%. 
The best fit point has Am 2 - 7.1 x 10~ 5 eV 2 and 
sin 20 — 0.83. Using the best fit point in Eq.lJH)! we 



obtain L P 



0.5 x 10- 5 T, 



-l 

MeV^MeV 



E 



The condition in 



Ea. (frp|) gives TueV < 2.8 and < 2.1 respectively for 
EMeV = 5 and 20. Similarly the resonance length L res ~ 
211 kilometers and 845 kilometers are obtained respec- 
tively for Ej^eV — 5 and 20. Using the resonance value of 
L e in Ea.lJT2|) we obtain M b ~ 1.15 x lO^R 3 . M Q E M \ V 
which is independent of the fireball temperature. For 
E M eV = 5 and 20 we obtain M h ~ 0.23 x 10~ 9 ] " 



^M 



and Mb ~ 0.58 x 10~ "R^Mq respectively. One can see 
that, in this case the temperature of the fireball is less 
compared to the lower limit of 3MeV. Also the baryon 
load in the fireball is less, which can be improved by 
increasing R 7 . By considering the solar neutrino mix- 
ing angle and the Am 2 , for neutrino oscillation in the 
fireball, the fireball temperature has to be less than 3 
MeV. If the L e vary with distance then at the resonance 
K r ~ 0.8 x 10~ 14 Z cm (£/ which implies, L' e — dL e /dx 
has to be very small to satisfy the condition in Ea. i|ll|) . 

The Super-Kamiokande Collaboration recently 
reported the atmospheric neutrino oscillation 
parameters 25] 1.9 x 10~ 3 eU 2 < Am 2 < 3.0 x 10" 3 eU 2 
and 0.9 < sin 2 29 < 1.0 with a 90% confidence level. 
Taking Am 2 



2.5 x 10~ 3 eU 2 and sin 2 29 - 0.9 we 



■1-3 

1 MeV ^ MeV 



e: 



and this implies 



obtain, L e ~ 0.98 x lO" 4 ^ 
T^MeV < 5 and < 3.8 respectively for EueV = 5 and 
20. So these temperature come within the range of the 
fireball as discussed earlier. The resonance length for 
EMeV = 5 an d 20 are respectively 5.2 kilometers and 21 
kilometers. Also for the above ranges of neutrino energy, 
the baryon load in the fireball is M b ~ 0.44 x 10~ 8 -RfM Q 
and M b ~ 0.11 x 10~ 8 i? 3 ,Af Q respectively. Thus before 
coming out of the fireball, the neutrinos can have many 
resonant oscillations. If the lepton asymmetry vary 
with distance, then l cm /L' e > 10 12 to have resonant 
conversion. 

Thirdly let us discuss the implication of the LSND 
and KARMEN results |26j on the neutrino oscillation in 
the GRB fireball. The combined analysis of LSND and 
KARMEN 2 results give 0.45 eU 2 < Am 2 < 1 eU 2 and 
2 x 10~ 3 < sin 2 29 < 7 x 10~ 3 with a confidence level of 
90%. We consider Am 2 - 0.5 eU 2 and sin 29 - 0.07 to 
estimate the fireball parameters. This gives, TMeV < 18 
and < 14 respectively for E^ieV — 5 and 20 and for 
these two values of EueV the baryon load is respectively 
M b - 0.3 x 10- 5 R^M Q and M b - 0.7 x IO^R^Mq. The 
resonance lengths for these two neutrino energies are re- 
spectively 0.4 kilometer and 1.4 kilometer. So this shows 
that, the propagating neutrinos have to oscillate several 
times, before they come out of the fireball. If the lepton 
asymmetry has a profile then in this case l cm /L' e > 10 10 
to satisfy the resonant oscillation. 

If the neutrino oscillation parameters are in the at- 
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mospheric and/or LSND range, the oscillation length at 
resonance within the fireball will vary from few meters 
to ~ 21 kilometers. As the resonance length is small 
compared to the size of the fireball, there will be many 
resonant oscillations for v e <-» before they emerge out 
of the fireball. So the average conversation probability in 
this case is V{t) <~ 0.5, where we have considered the fact 
that L e does not vary with distance. But if the oscilla- 
tion parameters are in the solar neutrino range, probably 
it is difficult to have neutrino oscillation. This is due to 
the fact that, the fireball temperature in this case is less 
(i.e. < 3 MeV). 

For the oscillation v e «-> v s , with charge neutral fireball 
and L p — L n , the effective potential is given by 



V ~ 4.02 x 10 



-12rp3 

1 MeV 



Le 

2 



-6.14x 10" y T 



MeV 



MeV. 
(17) 



Resonance condition can also be satisfied for the above 
process for V > 0. Both the processes v e <-> v^. T and 
v e ^> v s are equally probable if we consider the atmo- 
spheric neutrino oscillation parameter or the LSND one, 
with average probability ~ 0.5. Finally let us consider 



the process i/„ 



for which 



V 



2.0 x 10 _12 T; 



3 T 



MeV. 



(18) 



Here the higher order contribution to the neutrino poten- 
tial is absent and due to this, there is no constraint on 
the baryon loading of the fireball. Only the charge neu- 
trality of the fireball is sufficient enough for the resonant 
oscillation to take place. 



From the collapsar or merger model of SN with SN 
and/or from the collision of SN with a black-hole, lots 
of neutrinos will be produced and fractions of these neu- 
trinos will propagate through the fireball. Here we have 
studied the propagation of these neutrinos through the 
GRB fireball by assuming the later to be spherical with a 
radius R, charge neutral and L p = L n . Also we have as- 
sumed that L e > 6.14 x 10~ 9 T MeV so that the potential 
difference will be positive for neutrinos and there can be 
resonant oscillation. By using the known neutrino mass 
square difference and mixing angle from the solar, atmo- 
spheric and reactor experiments we estimate the fireball 
temperature, the baryon load and the lepton asymmetry 
in it. Our result shows that, if the neutrino oscillation pa- 
rameters are in the solar neutrino range, probably, there 
will be very few or no oscillation will take place. On the 
other hand if the neutrino oscillation parameters are in 
the atmospheric and/or LSND range as discussed above, 
there can be many resonant oscillation for v e <-> i^, r , 
v e <-> v s and v^^ <-> v a before they emerge out of the 
fireball and about 50% of these neutrinos will resonantly 
convert. These MeV neutrinos signals will be similar to 
the one from Supernovae, for example SN1987A. Unfor- 
tunately with the present generation neutrino telescopes 
these ^s can't be detected due to their cosmological dis- 
tance and the fluxes are extremely negligible compared 
to the galactic supernovae. 
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